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ABSTRACT 

We report polarization measurements in two prompt emissions of gamma-ray bursts, GRB 110301 A and 
GRB 1 10721 A, observed with the Gamma-ray burst polarimeter (GAP) aboard IKAROS solar sail mission. We 
detected linear polarization signals from each burst with polarization degree of II = 70 ± 22 % with statistical 
significance of 3.7 a for GRB 1 10301A, and II = 84^j % with 3.3 a confidence level for GRB 1 10721 A. We 
did not detect any significant change of polarization angle. These two events had shorter durations and dimmer 
bright ness compared with GRB 100826 A, which showed a significant change of polarization angle, as reported 
in lYonetoku et alj ( 1201 lb). Synchrotron emission model can be consistent with all the data of the three GRBs, 
while photospheric quasi-thermal emission model is not favorable. We suggest that magnetic field structures in 
the emission region are globally-ordered fields advected from the central engine. 

Subject headings: gamma-ray burst: individual (GRB 11 0301 A, GRB 1 10721 A) - instrumentation: polarime- 
ters - polarization - radiation mechanisms: non-thermal - relativistic processes 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) are the most energetic explo- 
sions in the universe. Since the discovery of GRBs in 
late 1960s, spacial distribution (directions and redshifts), 
lightcurves and energy spectra have been observed for large 
amount of GRBs. However, the emission mechanism of 
prompt GRBs, how to release the huge energy in short time 
duration, is still a crucial issue. The other characteristic of 
electro-magnetic wave, i.e. polarization in gamma-ray band, 
is thought to be a key to solve the problem. Polarization mea- 
surement can probe the existence of magnetic fields, and/or 
some kinds of geometrical structures, which are difficult to be 
examined from time variability and energy spectra. 

Several previo us works reported margin al detections of lin- 
ear polarization. Coburn & Boggs (2003) reported the detec- 
tion of strong polarization from GRB 021206 by RHESSI. 
However, independent authors analyzed the same RHESSI 
data, and con cluded that any polar i zation signals wer e 
not confirmed dRutledge & Foxl I2004t iWigger et alj 120041) 
KalemcietalJ 00071) : IMcGlvnn et alJ 120071) : IGotz etaD 
d2009j) reported marginal detections of polarization with ~ 2a 
confidence from GRB 041219 by INTEGRAL-SPl and -IBIS 
data. As the authors themselves discussed that the possibil- 
ity of instrumental systematics could not be completely re- 
moved, and a part of these results were obviously inconsis - 
tent with each other (see Section 1 of Yonetoku et al. 201 lb). 
Therefore, further observations with reliable instruments are 
strongly required to confirm the existence of polarization in 
prompt GRBs. 

Recently, Gamma-ray burst polarimeter (GAP; 
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aboard the IKAROS solar-power- sail foawaguc hi et alj 
2008; Mor i et all 120091) measured a linear polarization of 
n = 27 ± 1 1 % w ith 2.9 a confidence lev el in extremely bright 
GRB 100826A (lYonetoku et al.ll2011bl) . It was found that 
the polarization angle changes during the prompt emission 
with 3.5 a confidence level. The GAP detector is designed 
for the gamma-ray polarimetry of prompt GRBs, and also 
realized the quite small systematic uncertainty of 1.8 % level 
(lYonetoku et alJ [201 lal) . This may be the most convincing 
detection of polarization degree so far. But this is only one 
case, we need more sample to establish the existence of 
gamma-ray polarization in prompt GRBs. In this Letter, we 
report two more detections of gamma-ray polarization in 
GRB 11 0301 A and GRB 11 0721 A observed with IKAROS- 
GAP. We analyzed all data observed by GAP until the end of 
2011, and found significant polarization signals from these 
events. 

2. OBSERVATIONS 

We observed two GRBs, GRB 110301 A and 
GRB 110721A, with the Gamma-Ray Burst Polarimeter 
(GAP) aboard the IKAROS solar power sail. The polarization 
detection principle of GAP is to measure the anisotropic dis- 
tribution against the Compton scattering angle. According to 
the Klein-Nishina cross section, the gamma-ray photons tend 
to scatter toward the vertical direction of polarization vector. 
Therefore the modulation against the azimuth scattering 
angle is observed if the incident gamma-rays are polarized. 
The advantages of GAP are the high axial symmetry in shape 
and the high gain uniformity. These are the key to reduce the 
instrumental systema tic uncertaintie s . Deta ils of the GAP 
detector are shown in lYonetoku et al.l d201 lal) . 

The GAP detected GRB 11 0301 A on 2011 March 1 at 
05:05:34.9 (UT) at 0.946 AU away from the Earth. This burst 
was also detected by Fermi-GBM. The coordinate was deter- 
mined as (a, 8) = ( 229.35,-1-29 . 40) w ith an uncertainty of 1.0 
degrees in radius dFolev et al.ll201 ll) . which corresponds to 
48 ± 1 deg off-axis from the center of GAP field of view. Fig- 
ure Q] (top) shows the lightcurve of GRB 11 0301 A observed 
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with GAP in the energy range of 70-300 keV. 

Accor di ng to the time resolved spectral analyses by 
ILu et al.l d2012l) . the spectral peak energy (E p ) slightly 
changes during the burst, and show hard to soft trend from 
110 keV to 26 keV. Therefore GAP mainly observed the en- 
ergy range of E > E p . The energy fiuence in 10-1,000 keV is 
(3.65 ±0.03) x 10" 5 erg cm" 2 (Fol ey et alJ201 11) . 

GRB 1 10721 A was detected on 201 1 July 21 at 04:47:38.9 
(UT) at 0.699 AU from the Earth. Figure[l](bottom) shows the 
lightcurve of GRB 11072 1 A. This burst was first discovered 
by Fe rmi-GBM and LAT dTiernev et al.ll201 It IVasileiou et alj 
120111) . After that, the Swift-XRT pe rformed the follow- 
up observation of its X-ray afterglow dGreiner et alj 120 l"Tt 
Gruoe et alJl20Tlh . The coordinate is precisely measured as 
(a, S) = (333.66,-38.59) which corresponds to 30 deg off- 
axis. The optical t ransient was also detected by GROND 
dGreiner et al.l 1201 ll) and its redshift was measured as z = 
0.382 from tw o absorption lines of Call with Gemini-South 
dBerger et alJl20TTb . 

The spectral parameters, especially the E D value s dra- 
matically change dur i ng the burst dTiernev et alj 1201 It 
iGolenetskii et all 1201 It ILu et al.ll20T2h . The E p around the 
maximum intensity is about E p = 1130^|q keV, and the one 
of time integrated spectrum is E p = 393+[q4 keV. GAP mainly 
observed in the energy range of E < E p . The energy fiuence in 
10-1, 000 keV is (3.52 ± 0.03) x 10" 5 erg cm" 2 dTiernev et al.1 
120111) , which is very similar to GRB 1 10301 A. 

3. DATA ANALYSES 

3.1. Average Properties of Polarization 

We analyzed polarization data during the time intervals be- 
tween two dashed lines shown in Figure [T]for GRB 1 10301 A 
and GRB 11 0721 A, respectively. GAP obtained the polar- 
ization data between -16 s to 176 s since the GRB trigger. 
Since the time durations of these GRBs are relatively short, 
we used the background obtained in the same data. The net 
background rate for the polarization data is 60.0 counts s" 1 
for GRB 110301 A and 51.6 counts s" 1 for GRB 110721A. 
The total numbers of gamma-ray photons after subtracting 
the background are 1,820 and 1,092 photons for each burst, 
respectively. 

To estimate the systematic uncertainty, we first consider the 
spin rate of IKAROS spacecraft. The rotation of instrument 
generally reduces the systematic uncertainty because the dif- 
ferences of each sensor and the geometrical skewness are av- 
eraged. However, in these case, the time durations of bursts 
are smaller than the period of rotation of IKAROS space- 
craft. The spin rate is 1.61 rpm and 0.22 rpm for the epoch 
of GRB 11 0301 A and GRB 11 0721 A, respectively. Using 
the background interval of the data, we created the history 
of background modulation curves with the same time interval 
we analyzed, and confirmed each modulation was consistent 
with constant within the statistical error. We confirmed the 
systematic error due to the data analysis of short time dura- 
tion is about <7 JVJi i = 1 .0 % of the total polarization signals for 
each bin of both GRBs. 

Next, we estimated the systematic uncertainty between the 
detector response calculated by the Geant 4 simulator and the 
experimental data, which is mainly due to the off-axis direc- 
tion of incident gamma-ra ys. We performed seve ral ground 
experiments described in Yoneto ku et al.l (1201 lbl) with the 
proto-flight model of GAP. We estimated the systematic un- 
certainty was a sys 2 = 1 .9 % of the total polarization signals 



for each bin. 

In Figure [2] we show the modulation curve (polarization 
signals) after subtraction of the background. The error bars 
accompanying with the data (filled circles) includes not only 
the statistical error (a srar ) but also the systematic uncertainties 
described above. The total errors are calculated as af olal = 

°frar + a lys l + °fys 2 f° r eacn ^ m °f polarization data. 

The model functions (detector responses) were calculated 
with the Geant 4 simulato r consi dering the spectral evolu- 
tions reported by ILu et al.1 (120121) . who performed spectral 
analyses for 20 and 14 time intervals of GRB 110301 A and 
GRB 1 10721 A, respectively. Using their spectral parameters, 
we simulated the model functions for each time interval, and 
also combined into the one with the appropriate weighting 
factor estimated with the brightness histories. 

In these analyses, the free parameters are the polarization 
degrees (II) and angles (<p p ). We simulated the model func- 
tion with step resolutions of 5 % for polarization degree and 
5 deg for phase angles. In Figure [2] we show the best fit 
model with solid black lines, and also superposed the non- 
polarization model as the dashed lines on the same panel for 
easy comparison. The best fit parameters are II = 70 ± 22 % 
and 4> p = 73 ± 11 deg with y 2 = 14.0 for 10 degree of free- 
dom (d.o.f) for GRB 110301 A, and II = 84^j % and c/> p = 
160 ± 11 deg with X Z = 7.3 for 10 d.o.f for GRB 110721A, 
respectively. Here the quoted error are at 1 a confidence for 
two parameters of interest (n and <p p ), and the <fi p is measured 
counterclockwise from the celestial north. 

We show Ay 2 maps in the (n, cf> p ) plane in Figure [3] The 
white dots are the best fit results, and we calculate Ay 2 val- 
ues relative to these points. The la, 2a and 3a confidence 
contours for two parameters of interest are shown in the same 
figures. The null hypothesis (zero polarization degree) can be 
ruled out with 3.7cr (99.98 %) for GRB 110301 A, and 3.3a 
(99.91 %) confidence levels, respectively. Although these re- 
sults have relatively large errors compared with the previous 
GAP result for GRB 100826A (n = 27 ± 1 1 %, 2.9ct signifi- 
cance level), the polarization degree of these two GRBs may 
be larger than that of GRB 100826 A. From these observa- 
tions, we conclude the gamma-ray polarization really exists 
in the prompt GRBs. 

3.2. Time Resolved Polarization Analyses 

I Yonetoku et al.l (1201 lbl) found a significant change of po- 
larization angle in very bright GRB 100826 A. Therefore we 
also tried to investigate the existence of such effect in these 
two GRBs. 

As shown in Figure Q] (top), two main peaks exist in the 
lightcurve of GRB 1 10301 A. We divided the polarization data 
into two parts in time, i.e. s < t < 3 s and 3 s < t < 7 s. Here 
t is the time since trigger. The fitting results are n = 65 ±28 % 
and (j> p = 73 ± 14 deg for the first half part, and n = 80!^ % 
and <j) p = 7 1 ± 1 5 deg for the second half part. The significance 
of polarization detection is 2.87 a and 2.56 a, respectively. 

The brightest part of GRB 1 1 072 1 A lightcurve shows rather 
flatter compared with the standard "fast-rise and exponen- 
tial decay" of GRBs, and we divided the polarization data 
into s < t < 2 s and 2 s < f < lis. The best fit results 
are n = 8I+40 % and <\> p = 155 ± 14 deg for the first part and 
n = 78!^ % and cf> p = 164 ± 17 deg for the second part. The 
significance is 2.55 a and 2.16 a, respectively. 

We analyzed the polarization data of several time inter- 
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vals, but we could not find any significant change of polar- 
ization angle and degree. In these two cases, we conclude 
the polarization angles are stable during the prompt emis- 
sion within the err or. This is quite differen t properties from 
GRB 100826 A by lYonetoku et all (1201 lbh . and we should 
consider an explainable conditions of emission mechanism in- 
cluding both polarization properties. 

4. DISCUSSION 

Major results of the GAP observations so far are simply 
summarized as follows: (1) there are cases with and with- 
out a significant change of polarization angle (PA), and (2) 
GRB 100826 A, with a long duration T ~ 100 s, has a PA 
change, and its polarization degree (PD) is II = 27 ± 1 1 %, 
while GRB 1 10301 A and GRB 1 10721 A, with short durations 
T ~ 10 s, have no PA change, and their PDs are II > 30% at 2a 
error region (see Figure[3]l. Here we present some theoretical 
implications from these observational results for synchrotron 
and photospheric emission models, following Yonetoku et al. 
(I2011bh . 

First we discuss synchrotron models with different types of 
magnetic field structure. One type is helical magnetic fields, 
which may be advected from the central engine through the 
jet. S uch globally-ordered fields can pro duce high-PD emis- 
sion (iGranoj 120031: iLyutikov et all 120031) . We call this 'SO 
model' . In this model, the existence of a significant PA change 
implies that emission region consists of multiple patches with 
characteristic angular si ze D much smaller tha n the open- 
ing angle of the jet Oj (Yon etoku et al.l 1201 lbl) . Relativis- 
tic beaming effect only allows us to observe angular size of 
~ r -1 around the line of sight (LOS), where T is the bulk 
Lorentz factor of the jet. In the case of T" 1 ~ Oj, it is natu- 
ral that we see multiple patches with different magnetic field 
directions, and then we observe significant PA changes. On 
the other hand, if T" 1 <C Oj, we only see a limited range of 
the curved magnetic fields, which leads to no significant PA 
change even for patchy emission. The net PD of the latter 
case can be estimated as II ~ 40% for photon index Tb ~ 
-1 (whereas synchrotron emission from a point source has 
ri max = (-Fb)/(-Fb + §)) (see the case of yj = 100 in Figure 2 
of lToma et al. 2009), which is compatible to the observed PD 
of II > 30% from GRB 1 10301 A and GRB 1 10721 A. There- 
fore, the case of T -1 ~ 0j could correspond to GRB 100826 A, 
while the case of T~ l <C Oj to the other two bursts with no PA 
change. 

We may consider an alternative model in which GRB jets 
consist of multiple impulsive shells which have globally- 
ordered transverse magnetic fields with different direction for 
each shell. Let us call this 'non-steady SO model'. It has 
been recently claimed that suc h impulsive s hells can be ac- 
celerated to relativistic speeds dGranot et al.1120121) . This po- 
larization model may also apply to a scenario in which initial 
globally-ordered helical fields get distorted, making different 
field directions with in the angular size o f T" 1 , as considered 
in ICMART model dZhang & Yanll201 ll) . In this model, PA 
changes can naturally occur, although it does not necessarily 
occur when the number of shells is small so that T is rela- 
tively small. The net PD will be II ~ iEx/Vw, where N is 
the number of shells with different field directions. This could 
be as high as Il^ax ~ 60% for ~ -1 when T is relatively 
small, which is consistent with the observed PD of II > 30%. 

Shocks formed in the jet may produce sizable magnetic 
fields with random directions on plasma skin depth scales. 



Synchrotron emission from such fields can have high PD, pro- 
vided that the field directions are not isotropically random but 
random mainly in the plane parallel to the shock plan e (i.e., 
perpe ndicular to the local direction of expansion) dGranotl 
120031: iNakar et al.ll2003h . We call this 'SR model'. The linear 
polarization directions in the observer frame are symmetric 
around the LOS. Thus, if the emission i s patchy, net PD can 
be nonz ero and can have PA change s (lYonetoku et al.ll201 lbl) 
(see also Lazzati & Begelman 2009). The characteristic angu- 
lar size of the patches may be hydrodynamically constrained 
to be Op > T" 1 . In this model, the PD of emission from one 
patch strongly depends on the viewing angle of the patch V , 
i.e., the angle between the axis of the patch and the LOS, and 
it is limited as II < 40% (see the cases of y y > 1 in Figure 3 of 
iToma et al.ll2009l) . This implies that fine tuning, V ~ p + T~\ 
is required to have II > 30%. Furthermore, the bursts we ob- 
served are all very bright, which suggests that some patches 
are seen with V < P , for which II <C 40%. This suggests that 
the SR model is not favorable to explain the observed PD of 

n > 30%. 

Internal shocks may also produce strong magnetic fields 
with random directions on hydrod ynamic scales (llnoue et al.l 
l20Tlt iGruzinov & Waxmanlll999l) . We call this 'SH model', 
which produces net PD as II ~ ^lm"x/VN, where N is the 
number of independent patches with coherent magnetic field 
in the observable angular size ~ T" 1 , and can naturally lead 
to PA changes. Unlike the SR models, the emission from 
patches seen with small V has high PD, so that this model 
is in agreement with the high brightness of the bursts. How- 
ever, recent MHD simulations of i nternal shocks wit h initial 
density fluctuations imply N ~ 10 3 (Inoue et al. 201 1), which 
cannot explain the observed PD of n > 30%. 

Lastly, we discuss photospheric emission (Ph) model. This 
model assumes that the emission at E > E p is the quasi- 
thermal radiation from the photosphere. The emission at 
E < E p might be a superposition of many o f the quasi-thermal 
components wit h different temperatures dRvde et al.l 120101 : 
IToma et al]|201 ll) or contribution from synchrotron emission 
( Vurm et a l.l201 lb . The quasi-thermal radiation can have high 
PD when the radiation energy i s smaller than the baryon ki- 
netic energy at the photosphere (Beloborodo vl201 lh . The lin- 
ear polarization directions in the observer frame are symmet- 
ric around the LOS, same as the SR model. The PD of emis- 
sion from a given point may be determined by the brightest 
emission, coming from just below the photosphere, which can 
have II < nSax ~ 40% (lBeloborodovl201 ll) . The PD of emis- 
sion from one patch, whose angular size is hydrodynamically 
constrained to be p > T -1 , will reduce to n < 30% (in the 
same way as in the SR model). The patches with small view- 
ing angle V have II <C 30%. Therefore, the Ph model requires 
very fine tuning of V to reproduce the PD of GRB 1 10301 A, 
mainly at E > E p (i.e., dominated by the quasi-thermal com- 
ponent). 

To summarize, (1) the SR, SH and Ph models are not favor- 
able to reproduce all of our polarimetric observation results 
of the three GRBs, and (2) the SO and non-steady SO could 
explain all of them. 

Recently, early optical polarization from the forward shock 
of a GRB afterglow has been detected asll~ 10.4 ±2.5 % 
(at t = 149-706 s) for GRB 091208B (lUehara et al.ll20i2h . If 
the magnetic field directions in the emission region are ran- 
dom on plasma skin depth scales (i.e., the SR model for a 
shock produced in the circumburst medium), the observed po- 
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larization r eaches a max imum value around the jet break time 
(~ 1 day) ([Lazzati 2006|). The measurement of the early polar- 
ization, higher than the typ ical observed late-ti me polarization 
of ~ 1 -3 % (at t ~ 1 day: ICovino et al.ll2004h . thus disfavors 
this model. This result is consistent with our argument against 
the SR model for a shock produced within the jet. 

The Fermi satellite team has suggested that GRB 1 10721 A 
has a blackbod y component with te mperature ranging in the 
~ 10- 100 keV ( Axe lsson et al.l2012l) . which almost coincide 
with the GAP range 70-300 ke V. However, their fitting model 
of the integrated spectrum (their Figure 2) includes the black- 
body flux as only ~ 1 /5th of the total flux in the GAP range. 
Thus reduction of PD by addition of the blackbody compo- 
nent is small. The polarization degree and angle in the GAP 
range is practically determined by those of the non-thermal 



component, which could be synchrotron emission. 

The SO models assume globally-ordered fields in the emis- 
sion region. On the other hand, observations have suggested 
that prompt emission has ver y high efficiency (even > 90% 
for some bursts: IZhang et all 120071: lloka et all 120061) . which 
means that energy dissipation, usually involving field distor- 
tion, occur globally. Reconciling high PD with high efficiency 
looks a dilemma, which will have to be resolved in more quan- 
titative modeling. 
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FIG. 1. — Lightcurves of the prompt gamma-ray emission of GRB 110301A (top) and GRB 110721A (bottom) detected by GAR The vertical dashed lines 
indicate the time interval of polarization analyses for each burst. 
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FIG. 2. — Number of coincidence gamma-ray photons (polarization signals) against the scattering angle of GRB 1 10301A (top) and GRB 1 10721A (bottom) 
measured by the GAP in 70-300 keV band. Black filled circles are the angular distributions of Compton scattered gamma-rays after the background subtraction. 
The solid lines are the best fit models of each event calculated with the Geant4 Monte-Carlo simulations. The model functions of non-polarized cases are 
superposed on each panel for easy comparison. 
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FIG. 3. — Two A\ 2 maps of confidence contours in the (JL,(j>p) plane for GRB 110301A (top) and GRB 110721A (bottom). The white dots are the best-fit 
parameters, and we calculate A\ 2 values relative to this point. A color scale bar along the right side of the contour shows levels of A\ 2 values. The significances 
of the existence of linear polarization are 3.7 a for GRB 1 10301A and 3.3 cr for GRB1 10721A, respectively. 



